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Mitgliedsbeitrag fur 1992 


Die DUG-Mitglieder aus Osterreich werden gebeten, den Mitglieds- 
beitrag von S 300.- fur das Jahr 1992 mit beiliegendem Zahlschein zu 
begleichen ! 

Die Mitglieder von auLerhalb Osterreichs bitten wir ihren Beitrag 
fur 1992 uber DM 45.- nach Moglichkeit in Form eines Euroschecks an 
den Herausgeber zu schicken. Alle anderen Uberweisungsf ormen verur- 
sachen zum Teil recht hohe (bis zu DM 20.-) Bankspesen. 

Herzlichen Dank 


Membership Dues for 1992 


The Austrian DUG Members are asked for settling the memership dues for 1992 using the enclosed 
paying-in form! 

We invite the members from outside Austria to pay their dues of GM 45.- for 1992 using an euro- 
cheque as far as possible and sending it to the editor. Each other kind of payment causes considerable 
banking charges (sometimes GM 20.-). 


Thank you 





D-N-L#4 


ir t t r u or THr rviron 


p i 


Lieber Derive Anwender! 

Mit der 4. Ausgabe des D-N-L beendet die 
DUG das erste Jahr ihres Bestehens. Die wei- 
terhin steigende Mitgliederzahl - derzeit 205 - 
ist ebenso erfreulich wie die Tatsache, dass uns 
immer wieder interessante Artikel erreichen. 
Trotzdem bitten wir um weitere Beitrage aus 
alien Gebieten der Lehre und Anwendung. Ich 
mochte besonders die technischen Anwender 
auffordern, uns iiber ihre Einsatzmoglichkeiten 
von DERIVE zu berichten. Schreiben Sie auch, 
wenn Sie Ideen haben, die Sie geme be- 
handelt wissen oder zur Diskussion gestellt 
haben wollen. 

Auf Seite 37 finden Sie die Ausschreibung 
fur eine D£/?/K£-Fruhjahrs-Schule in Krems, 
einer netten Stadt in Niederosterreich. Wie Sie 
vielleicht wissen, hat das dafiir zustandige 
Ministerium fur einen GroBteil der Hoheren 
Schulen eine Generallizenz fiir DERIVE er- 
worben. Die damit befassten Herren und da 
ganz besonders LSI Dr. Heugl und MR Dr. 
Szirucsek erkennen die Herausforderung, die 
die Computeralgebra fiir die Schulmathematik 
darstellt und wollen ihr offensiv begegnen. 
Man darf entscheidende Impulse von dieser 
Tagung erwarten. 

Wie schon im letzten D-N-L angekiindigt, 
iiberschreitet die DUG mit dem Jahr 1992 
offiziell die Grenzen Europas und ladt auch die 
US-DERIVE-Usqt zum Beitritt ein. Wir kon- 
nen uns schon jetzt auf den Ozean iiberspan- 
nende Kontakte freuen. 

Zum SchluB will ich Ihnen noch ein frohliches 
Weihnachtsfest sowie ein gliickliches Jahr 
1992 in Gesundheit und Frieden wiinschen. 
Nehmen Sie die beigelegte Diskette mit alien 
.MTH-Dateien aus D-N-L# 1 bis D-N-L#4 als 
kleine Aufmerksamkeit fiir Ihr Interesse an der 
DUG, als Dank fiir Ihre Mitarbeit und zugleich 
als Bitte, auch weiterhin am Ausbau der DUG 
mitzuwirken. 


Dear Derive User, 

the DUG ends its first year of existence with 
the first issue of its bulletin, the D-N-L. The 
steadily increasing number of members is as 
pleasing —205 till now - as the fact that inter- 
esting contributions are reaching the editor's 
desk again and again. Nevertheless we ask for 
further papers from all fields of teaching and 
application mathematics. I'd like to invite es- 
pecially the technical users to inform us about 
their use of DERIVE. Please also write if you 
have any ideas which you want to be dealt with 
or you want to have discussed. 

On page 37 you will find the announcement 
for a DERIVE Spring School in Krems, a nice 
town in Lower Aus-tria. Perhaps you know, 
that the Ministry of Education has bought the 
general DERIVE licence for a great part of 
Austria's Secondary Schools. The people being 
occupied with this task and especially LSI Dr. 
Heugl and MR Dr. Szirucsek recognize the big 
challenge which Computer Algebra performs 
for maths in school and they will face this 
challenge bold-ly. We may expect a lot of im- 
portant impulses from that meeting. 

As announced in the last D-N-L, the DUG 
will go beyond Europe's borders in 1992 and 
will invite the US DERIVE Users for member- 
ship. We are looking forward to having con- 
tacts reaching over the ocean in the future. 

At last I wish you a Merry Christmas and a 
Happy New Year 1992 in health and peace. 
Take the disk enclosed containing all the 
.MTH-files from D-N-L# 1 to D-N-L#4 as a 
little present for your interest and cooperation 
and at the same time as an invitation to partici- 
pate in the further development of the DUG 

With my best regards 


Mit den besten GriiBen 
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The Derive-News-Letter is the 
Bulletin of the Derive-User- 
Group. It is published at least 
three times a year with a con- 
tents of 30 pages minimum. 

The goals of the D-N-L are to 
enable the exchange of experi- 
ences made with Derive as well 
as to create a group to discuss 
the possibilities of new me- 
thodical and didactic manners 
in teaching Mathema-tics . 

Editor : 

Mag. Josef Bohm 
A-3042 Wurmla 
D ’ Lust 1 

Phone 43- (0) 2275/8207 
Austria 


Contributions : 

Please send all contributions 
to the above address. Non-Eng- 
lish speakers are encouraged to 
write their contributions in 
English to reinforce the inter- 
national touch of D-N-L. It 
must be said, though, that non- 
English articles will be warmly 
welcomed none-the-less . Your 
contributions will be edited 
but not assessed. By submitting 
articles the author gives his 
consent for reprin-ting it in 
D-N-L. The more contributions 
you will send to the Editor, 
the more lively and richer in 
contents the Derive-News-Letter 
will be. 


Preview: Contributions for the next issues: 

Finding a gradient; 

Zur Problematik des Einsatzes von DERIVE; 
Cyclomania (Hypotrochoids & Epitrochoids); 

Two Underground Tunnels (applied Geometry); 
Module on sequences and series with Derive (Part 2) 


(will be published March 92) 


Impressum : 

Medieninhaber : DERLVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrif t 

Herausgeber: Mag. Josef Bohm, Adresse wie oben 
Herstellung: eigene Vervielf altigung 
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Waclaw Sadowski, Bukowiec 

Dear Sir, 

My name is Waclaw Sadowski and i am thirty.three years of age. I live in a small village called Male Konskie in 
the central region of Poland. I am a higher school graduate and I am a lecturer of mathematics in Secondary 
School. 

I am very much interested in getting know more on the Derive User Group and its bulletin and that’s why I’ve 
made up my mind to write to you. Recently I have been successful in entering into contacts with some scientists 
who involve in research on new developments in mathematics and mathematics teaching. 

The secondary schools in Poland are secondary schools which provide comprehensive secondary education. The 
graduates of this school, after having finished the last (fourth) grade, obtain certificates of secondary education 
and after having passed the final examination (so-called ’maturity’) are entitled to apply for admission to the 
institutions of higher education. The four-year secondary general schools realize plans and curricula of teaching 
approved by the Ministry of National Education. 

We Polish scientists need help from specialists from western countries to help Poland in its new form to integrate 
into the western part of Europe. I shall much appreciate if you could kindly take action to provide me with in- 
formation on Derive and one free of charge your bulletin. I look forward to hearing from you before long. 

Yours sincerely 


Madia Lokar, Kranj 

Dear Editor! 

Few days ago I got the Derive news letter# 1. 1 found it very interesting and I would like to join the Derive user 
group. I would do that immediately as I manage the way to pay Annual Membership Dues, as it is now almost 
impossible to do, because of the situation in our country. 

[..some comments to the equation in D-N-L#l, which leads to x7 - 1 = 0] 

Although I mostly work with Mathematica, I am very interested in Derive and found it very good. In the last 
half a year, me and a couple of friends tried to put Derive on the list of software, specially convenient for use in 
secondary schools. We hope we manage that and to the end of year it is expected that all of 135 secondary 
schools will receive a copy of Derive for inspection. We also plan to establish special courses about Derive. 
Therefor I ask you for possible help. Any suggestions, experience with training for use of Derive would be appe- 
ciated. 


With best regards 


D-N-L: We are very glad that the DUG might reach the young democratic states in Eastern Europe. If there are 
any DUG members who want to establish contact with Mr. Sadowski or Mr./Mrs. Lokar, then please write or 
call for their addresses. 


Herbert Appel, Schweinfurt 

[...] Zusatzlich erhalten Sie von mir mit der Bitte um Veroffentlichung einen weiteren Artikel zur Problematik 
des Einsatzes von DERIVE im Mathematikunterricht. Der Artikel kann, falls notwendig, sinnvoll gekurzt wer- 
den; ein vollstandiger Abdruck (vielleicht schon im DNL#4) ist mir jedoch am liebsten. 
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Verbesserungsvorschlage zum DNL: 

- konnte dicker sein. 180 Mitglieder miiBten eigentlich eine ganze Menge Artikel produzieren. Ich ware unter 
Umstanden auch bereit, einen hoheren Mitgliedsbeitrag zu bezahlen. 

- Bitte an die Mitglieder, bei Einverstandnis der anderen Mitglieder, Namen und Adressen weitergeben, um 
miteinander in Kontakt treten zu konnen. 

- Ich benutze Derive auf meinem Atari ST mit ATonce 16Mhz. Derive spricht hier im Gegensatz zu "Cabri- 
Geometrie", die VGA-Graphikkarte nicht an. AuBerdem kann ich trotz GRAPHICS aus DOS 4.0 mit PRTSCR 
keine Hardcopy auf meinem NEC P6 machen. 

- Bitte nicht mit Klammem beim DNL sparen. DNL # 2 war noch mit 3 Klammern am Auseinanderfallen gehin- 
dert. (DNL#3 nur noch eine Klammer). 

With all - autre chauses - sono contento - weiter so - do swidanje 
Im Rahmen der Volkerverstandigung. 


D-N-L: What is the opinion of the DUG-members? In the next issue we will publish the complete list of the 
DUG, exept those who have informed me not to be mentioned by name. 

I hope that the 202 members of the DUG are as productive as Mr. Appel thinks them to be. 

The bad finishing of DNL # 3 was caused by lack of time. Sorry, hope to do better now! 

Mr. Appel has sent an article concerning the problems of using DERIVE when teaching mathematics. It will be 
published in DNL# 5. 


Robert Setif, Annemasse 

1 . Pourquoi le chronometrage de certains calculs dans un environnement pratiquement identique donne-t-il 
parfois des resultats tres differents? Voir la fonction mod. 

2. La fonction mod de DNL# 2, p 15 m'a surpris. Elle ne fonctionne pas avec par exemple 

vector (mod (k, 2) , k, 0, 24) Simplify or Approx 

3. Voudriez-vous des fonctions aussi simples que MOD IT (k,n) : ^iterate (a-n, a, k, k/n) ou des 
fonctions ADD MUL DIV recursives tres peu performantes, mais interessantes pour apprentis- 
programmeurs? 

4. Si cela peut etre utile a quelqu’un, je peux poursuivre le calcul (a la main!) des fractions continues au -dela 
de a/700, mais je prefererais qu'on me trouve une fonction de Derive (ou de MuMATH) qui trouve 
automatiquement les resultats. 

5. Pourriez-vous approfondir le sujet des erreurs pour les fractions continues successives? C’est ma partie 
faible ! 

6. Auriez-vous d'autres sujets a me proposer qui ne demandent pas un trop grand quotient intellectuel, car je 
commence a etre .. use..(et je travaille tout en gardant mon petit fils de 8 mois!)? 
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7. Je ne suis pas equipe pour le bulletin electronique for Derive users. Pourriez-vous publier dans DNL les 
reponses les plus interessantes ? 

8. Y a-t-il une nouvelle version de Derive en chantier ? Avec quelles ameliorations ? 

9. Etes-vous utilisateur (ou 1 avez vous ete) de Mumath ? Des fichiers Mumath vous interessent-ils ? 

A mon avis, Mumath est plus interessant que Derive pour les calculs sur N et Q. 


Souhaits pour nouvelles versions de Derive 

1. Option notation: Affichage, par tranches de 3 chiffres. 

Par exemple 2 8 * 10 = 1 048 576 

a/2 = 1.414 213 562... 

2. Option notation: affichage dans tout 1’. ecran. 

3. Option mute: Signal sonore a la fin d'un calcul, car c'est penible d’attendre (meme si pendant ce temps je 
peus calculer des fractions continues!) des resultats qui viennent apres plus de 1 heure ou meme plusieurs 
minutes. 

4. Option break: (comme dans MuMATH) qui permette d’interroger certaines variables, de comprendre des 
erreurs, de connaitre la situation quand on est "dans le noir". 

5. Plus de possibility de programmation , bien qu'il y ait d’ja un progres considerable depuis Derive 1., et ceci 
avec tres peu de mots nouveaux (Vector, Iterate, if). 

6. Passage Mumath <— > Derive comme il y a Fortran <— > Derive, au moins pour les resultats numeriques et 
les polynomes. 

7. Introduction de LISTES : Plus souples que les vecteurs, avec pop et push et liste de listes qui n’auraient pas 
forcement la meme dimension. 

8. Possibility d’une nouvelle sorte de fonction qui afficherait les resultats a mesure a chaque appui sur une 
touche speciale (par exemple "space" ASCII 32 ) au lieu d’attendre que par exemple le vecteur soit tout 
calcule avec les risques de "memory full". 

Avec mes meilleurs sentiments 


Bon vacances 

Merci pour les ameliorisations de mon fielder SYRACUSE. 





p 1 2 3 4 5 6 7 


Derive - User - Forum 


D-N-L#4 


Here you can find some improvements to the SYRACUS.MTH file from D-N-L#3. 


In 1991 it was necessary to produce one's own mod-function mod (x, y) , (see DNL#3). 


sy(n) := 

If M0D(n, 2) 

#1 : n/2 

3 ■ n + 1 

#2: auxsyr(xO) := ITERATESf IF(x > 1, sy(x), x) , x, x0) 

#3: auxsyr(3) 

#4: [3, 10, 5, 1 6 , 8, 4, 2, 1 , 1 ] 


|lsyr(n) := DIM(auxsyr(n) ) - 1 


#6 : lsyr(3) 

Length of the sequence 
# 7 : 8 

syr(x0) := UECTORC (auxsyr(xO) ) , i, lsyr(xO)) 

#8: i 

#9: syr(3) 

The complete sequence 

#10: [3, 10, 5, 16, 8, 4, 2, 1] 

A very long sequence with an ordinary start: 

#11: syr(27) = [27, 82, 41, 124, 62, 31, 94, 47, 142, 71, 214, 107, 322, 161, 484, 242, 121, 

364, 182, 91, 274, 137, 412, 206 , 103, 310, 155, 466, 233, 700, 350, 175, 526, 263, 790, 
395, 1186, 593, 1780, 890, 445, 1336, 668, 334, 167, 502, 251, 754, 377, 1132, 566, 283, 
850, 425, 1276, 638, 319, 958, 479, 1438, 719, 2158, 1079, 3238, 1619, 4858, 2429, 7288, 
3644, 1822, 911, 2734, 1367, 4102, 2051, 6154, 3077, 9232, 4616, 2308, 1154, 577, 1732, 
866, 433, 1300, 650, 325, 976, 488, 244, 122, 61, 184, 92, 46, 23, 70, 35, 106, 53, 160, 
80, 40, 20, 10, 5, 16, 8, 4, 2, 1] 


# 12 : 


maxsy(n) := ITERATES( 


r*v 


sy(u ), MAX( 


u , u )1 , u, [ 1 , n , 1 
2 3 J 


], lsyr(n) - 1) 


#13: maxsy(3) = 


1 3 1 

2 10 3 

3 5 10 

4 16 10 

5 8 16 

6 4 16 

7 2 16 


8 


1 


16 
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max_el(n) := (maxsy(n)) 

#14: lsyr(n),3 

#15: max_el(3) = 16 


Today it is easier than in 1991! 

#16: MAX(syr(3)) = 16 

#17: maxel(27) = 9232 

9232 is the maximum element of the sequence starting with 27 

#18: syracusfu, u) := TABLEf [lsyr(x) , max_el(x)], x, u, u) 


syracus( 1 , 

50) 

24 

11 

24 

25 

24 

38 

26 

11 

40 

27 

112 

9232 


#21: syrgr(u, u) := TABLE(max_el(x) , x, u, u) 

#22: syrgr( 1 , 50) 



LSYRGRCu, v) := VECTOR C[x , LSYRCx)], x, u, v) 
LSYRGRCL, 50) 

-120 



2 4 6 S 10 12 14 16 IS 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50 
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This time monsieur Setif opens his treasure box to show us some spirals and polygons. Until expression #5 work 
with rectangular coordinates, then switch to polar coordinates. (More graphs on p 38) 


# 1 : 
# 2 : 
#3: 
# 4 : 
#5: 
#6 : 
#7: 


SPIRR{r, k, a, 
SPIRR ^3, 0.95, 
SPIRR ^3, 0.95, 
SPIRR|l, 1.02, 
SPIRR|l, 1.05, 
SPIRR|l, 1.02, 
SPIRR 1, 1.02, 


a_, n) := ITERATES k-u 



0, , 50 

5 


2 n 

0, , 24 


2 n 

0, — - — , 30 
2 n 

0, , 45 




#12: POLVGON(r, a, n) := ITERATES u 


0 , 


2 n 


u, [r, a], n 


#13: P0LVG0N{1, 0, 5) 

#14: P0LVG0N{2, 0, 3) 

#15: STARS(r) := ITERATES u 


0 , 


4 n 

5 J' 


u, [r, 0], 5 


Try to find the N-Star! 


#16: STAR5{2) 

#17: SPIR(x0, y0, x_, v_, n) := ITERATES(u + [x_, y_] , u, [x0, y0], n) 

#18: SPIR{1, 0, 0.05, 90°, 32) 

f n 

#19: SPIR 3, 0, -0.05, , 60 

[ 2 

#20: SPIR{3, 0, -0.1, 120 °, 30) 

#21: SPIR{0.2, 0, 0.005, 30°, 288) 



SEGPOLREC(r, r_, a, a_, n) := 

If n = 0 

#26: [ [r, a]] 

APPEND{ [ [r + n r_, a + n a_] ] , SEGPOLREC(r, 
#27: SEGPOLREC 1, 0.1, 0, — , 15 

l 8 



r_, a, a_, n - 1)) 


#28: SEGPOLREC 0.5, 0.1, — , — — , 11 

l 4 12 J 

#29: U ECT OR { POL V GON { a , 0, 5), a, 1, 2.4, 0.2) 

#30: UECTOR(STAR5(a), a, 1.5, 2.5, 0.5) 
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Die Hilfsfunktionen konnen wir auch gleich in den ITERATES-Befehl mit aufnehmen: 
ITERATESCI" Fl(x , x ), F2(x , x ) 1 , x , fx , x 1 , n) 

L 12 1 2 J L 1 21 


Es ist gegebenenfalls zweckmaBig, sich dies nun in einem Hilfsfile ’’ITERATIO.MTH” abzuspeichern, 
das etwa folgendermaBen aussehen konnte: 

(In einer MTH-Datei konnen keine Textboxes gespeichert werden, Kommentare wie in #1, #2 usw. 
miissen unter " geschrieben werden.) 


It makes sense to save this in a utility file ITERATIO.MTH, which could look like as follows: 

(In a MTH-file textboxes cannot be saved, so include comments (#1 , #2, etc) under quotes ".) 

#1: file ITERATIO.MTH 

#2: The follwing system will be iterated: 

#3: xl_new=Fl(xl_old , x2_old) 

#4: x2_new=F2(xl_old , x2_old) 

#5: Define functions FI Cxi, x 2) and F2(xl,x2) 

#6: Call result_matr ix ([3 , 4] , 5) , with [3,4] is the initial vector x 

#7: and you want to get the first 5 iterations 

# S : ============================================================ 

#9: [FI (xl, x2) :=, F2(xl, x2) :=] 

result_matrix (xO , n) := ITERATES (T FI (x , x ), F2(x , x )1, x, xO, n) 

#10: L 1 2 1 2 J 

# 11 : ============================================================ 

#12: 1. Beispiel / 1. Example 

#13: [FI (xl, x2) := xl + x2, F2(xl, x2) := xl-x2] 


#14: 


result_matrix ([1 , 2], 3) 


' 1 2 
3 2 

5 6 

_ 11 30 _ 


3.3 Beispiel einer Anwendung 

Volterra-Zyklus zwischen Rauber- und Beutetieren Predator - Prey - Cycles 


Fur Rauber - Beute - Simulationen gilt das Gleichungssystem: 
Predator-Prey-simulations can be described by: 


jc/ =x l -(l-d + c-x 2 ) 
= x 2 -(l + a-b-x l ) 


mit xi = Rauberzahl und x 2 = Beutezahl (xi = predators, x 2 = prey) 


Dabei haben die Parameter die folgende Bedeutung: 

a Wachstumsrate der Beutetiere (growth rate of prey) 

b Verlustrate der Beutetiere (loss rate of prey) 

c Regenerationsrate der Rauber (rate of regeneration of predators) 

d Energieverlust der Rauber (loss rate of energy of predators) 
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Nachdem wir die Datei ITERATIO . MTH geladen haben, geben wir ein: 


#15: 

[a := 1 

, b := 0.002, c := 

0. 

00001 

#16: 

Fl(xl, 

x2) := xl ■ (1 - d 

+ 

0x2) 

#17: 

F2(xl , 

x2) := x2 ■ (1 + a 

- 

b-xl) 

#18: 

resultjnatrix ([350 , 6000] 

, 4) 


350 

6000 

343 

7800 

342.314 

10249.2 

350.013 

13481.5 

369 . 199 

17525.6 


#19: 


und erhalten nach SIMPLIFY die nebenstehende Matrix. 

3.4 Graphische Darstellung 

Fur unser Iterationsproblem haben wir in Abschnitt 3.3 die Losung in result_matrix (xO , n) 
zusammengefasst. Sie lieferte eine Matrix in der Form 


e : = 


" e u 

®12 


x i 

X 2 

®21 

e 22 


t 

x i 

! 

X 2 

e Bl 

e 32 


ft 

ft 

e 41 

e 42 

= 

x i 

X 2 




ttt 

ttt 




x i 

X 2 

_®n+l,l 

®n+l,2_ 





Wenn man die entwickelten einzelnen Komponenten von x\ und x 2 in Abhangigkeit von der Iterations- 
tiefe graphisch darstellen will, muss man daraus, wie in Abschnitt 2.4 besprochen, Matrizen in folgen- 
der Form erzeugen: 


pi c_l : = 


" 0 e il " 


" 0 e 12 

1 e 21 


1 e 22 

2 e 31 


2 e 32 

3 e 41 

und pi c_2 : = 

3 e 42 

n e . . 

_ n+l,l_ 


n e , _ 

_ n+l,2_ 


Mitm SUB i SUB j = - m ist die Matrix result_matrix (xO , n) - lassen sich beide 

Bilder mit dem VECTOR-Befehl durch eine Schleife einlesen. 


Expression to extract data for plotting the devel- 
opment of predator population (pic_1 ) and prey 
population (pic_2). 


p-ic_l := VECT0R(|” i - 1, m 




p-ic_2 := VECT0R(T i - 1, m 


Wir priifen dies nach, indem wir eine Matrix m eingeben: 


i,2 

1 2 
3 2 

5 6 

L 11 30 J 


i, DIM (m)) 
i, DIM (m)) 


3 11 


3 30 


... und erhalten richtig nach SIMPLIFY: 


[picJL, pic_2] = 
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Die Daten fiir unsere Simulation sind / Data for our simulation are 


[a := 1, b := 0.002,c := 0. 00001, d := 0.08,n := 100, xO := [450, 8050]] 


pi c_l(x0 , n) 


■1400 

■1200 

■1000 

■S00 

■600 



■400 

■200 

-200 


Die Rauber ■ 

The Predators " \ 


10 IS 20 25 30 35 40 45 50 55 60 65 70 75 SO S5 90 95 100 3 


c_2(x0,n) 


50000 

■ 

■ 

40000 

Die Beutepopulation 
The Prey population 

30000 

■ ■ 

20000 

■■ ■ 
■ ■ 

■ ■ 
■ 

■ 

■o 

■g 

5L 

■ 

. 

■ ..... 

■ ■ ■ 



5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 1 


Eine andere Wahl fiir die Parameter liefert eine mehr oder weniger zyklischen Verlauf der Populati- 
onsentwicklung / Another choice of parameters shows a cyclic up and down of both populations. 
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Anmerkung des Herausgebers: 

In der Zeitschrift "Spektrum der Wissenschaft" (Scientific American) fand ich in 2/84 und 2/85 jeweils 
einen Artikel uber das, von Herm Schumm bearbeitete Thema. Dabei wurden zwei interessante Gra- 
phiken gegenubergestellt: 



Furs of Snowhare and Lynx caught in Canada between 1850 and 1900 (above) 
and solution of the associated differential equation (below). 

Pictures from "Scientific American", 1984 and 1985. 

Dies ist die exakte Losung der zugehorigen Differentialgleichungen: 



Bild 5: Zyklische Schwankungen in der Population von Riiuber und Beute, 
die sich aus der Lotke-Volterra-Gleichung ergeben. 
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I enclose a short article on the use of DERIVE to find inverse Laplace Transform s which may be of use 
for the Bulletin of the Derive User Group. I have deliberately left it open-ended with some sugges- 
tions for further investigations. 

Yours faithfully 
J.V. Royle 

Inverse Laplace Transforms using DERIVE 

J.V. Royle, Leicester 


We wish to find the inverse Laplace Transform, L 


-i 


Pjs) 

kQ(s); 


P(s) 

of a rational function , where 

Q(s) 


deg Pis) < deg Q(s). The usual elementary method is to 


(i) 

(ii) 

(iii) 


factorise Q(s) into real linear and quadratic factors, 
P(s) 


express 


Q(s) 


in terms of Partial Fractions, 


use tables of standard Laplace Transforms to find the inverse Laplace Transform of each of 
the Partial Fractions. The Expand command in DERIVE can be used to find the Partial Fraction 
expansion required in step (ii), provided the denominator factorises into real linear and quad- 
ratic factors with rational coefficients. 


Example Find L 


f s 2 +4?-2 ^ 

y 2 +11/5' +1 8 s + 9 j 


s + 4-s - 2 


# 1 : 


# 2 : 


3 2 

2-s + 11 -s + lS-s + 9 


23 


3 ■ (2 ■ s + 3) 6 ■ (s + 3) 2-(s + 1) 


Simplif y>Expand delivers the factorisation. There is no need to first factorise the denominator. 

f 1 ^ 


Using the standard inverse Laplace Transform L 


\s — aj 


■ e at , we get 


s + As -2 


2V+11V +18^ + 9 


23 e 


5 5 e 


Among the User Contributed files you can find LaplaceTransforms.dfw which according to its header 
computes the inverse Laplace transform of the function f(s) in terms of t. The author of this file is our 
friend Terence Etchells, who revised it first time in August 2003. 

So we load this file as a Utility file and try to find the Inverse Laplace Transform of the above expres- 
sion to check this calculation - and to check the file, of course. 


In v Laplace 


+ 4* s - 2 


3 2 

2*5 + 11* s + 18* s + 9 


-t 


5*e 


23 *e 


- 3 * t/2 


3 * t 


5*e 
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Problems with this method arise if 

(i) there are Partial Fractions with quadratic denominators, 

2s + 5 

e.g. for — DERIVE gives the Partial Fraction expansion 

s + 6s + 16s + 16 


EXPAND 


4-(s + 


^ s + 6-s + 16-s + 16 


4 ■ (s + 4-s + S) 


2-fs + 4-s + SI 


Using standard inverse Laplace Transforms of the form 
( s + a } 

L 1 = e at cos (bt) still leaves some manipulation to be done which is not easily 

^ (s + cl ) + b j 

performed with DERIVE. 

(ii) the denominator, Q(s), does not factorise into linear and quadratic factors with rational coeffi- 
cients. 

These problems can be surmounted by the following method, which depends on the result if Q(s ) has a 

P(s) A 

simple factor of (s - a ) then the Partial Fraction expansion of contains a term , where 

Q(s) s-a 

A _m_ 

Q'(a ) ' 

Hence if Q(s ) is an n th degree polynomial with n simple roots a u a 2 , a m real or complex, and 
deg P(s) < n then: 


00 ) Tfs-a, 


where A. = 


Q\a t ) 


, then using L 


= e at the final result is 


P{s) 


Q(s)J t f&W 


7 - 


P(^) e a,t 


f P(s) | 

Hence to find L~ l , where Q(s ) has only simple roots, the procedure in DERIVE is: 

ISO) J 


(i) define a function ILT (a) : =LIM (p/dif (q, s ) A s , a) *exp (a • t ) 

(ii) author the numerator and denominator as variables p and q 

(iii) Solve q = 0, best to use Mixed Mode 

(iv) assuming there are as many different solutions, a u a 2 , a n as the degree of Q(s ) then the 

roots are all simple, find the inverse Laplace Transform by authoring and simplifying 


ILT(ai) + ILT (a 2 ) + ... + ILT(a n ). 
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P(s) 

There is no problem with complex roots since, if ^ J is real, complex roots occur in complex conju- 
gate pairs and DERIVE automatically expresses complex exponentials, which will occur in the ILT- 
sum as complex conjugate pairs, in terms of real trigonometric functions. 

In the following examples it is assumed that the ILT-function has been defined as above and that 

Precision Mode has been set to Mixed. 


Example 


Find L 


3s + 5 


s 3 + 6s 2 + 16s + 16 


#6: p := 3-s + 5 

3 2 

#7: q := s + 6-s + 16- s + 16 

3 2 

#S: s + 6-s + 16- s + 16 

#9: SOLVE (q, s) 

gives the result 

#10: [s = -2, s = -2 - 2 ■ i_ , s = -2 + 2-i_] 

so to find the inverse Laplace transform: 

#11: ILT (-2) + ILTC-2 - 2-0 + ILT(-2 +2-0 

gives the result 


2 ■ t 


f C0S(2 ■ t) 3 ■ SIN(2 ■ t) 


#12: e 


1 ^ 
4 j 


We check again applying InvLaplace: 

f 

InvLaplace 


' 3-5+5 

- 2 * t 

= e 

' C0S(2 • t) 3 ■ SIN(2 ■ t) 1 ' 

■ 

3 2 

1 

L 4 2 4 J 

+6-s + 16-s + 16 > 




So in fact this tool seems to be very reliable! 


Example Find L 


r s 2 -4s + 2 
s ^ + 3s 2 + 5s + 7 


#13: p := s - 4-S + 2 

3 2 

#14: q:zs + 3-s + 5-s + 7 

#15: SOLVE (q, s) 


#16: 


( 2 - J 57 . 2 1 

+ J 

p-* 7 - 2 ] 

l 9 J 

1 + 1 

l 9 J 

3-J87 7 ' 

r\i 

r 7 3 -js7 

S 2 J 

i + 1 

l 2 8 


C JS7 1 \L/3 

f JS7 

1 "jl/3 

l 36 

4 J 

l 56 

4 J 

r 787 1 y./3 1 

' 78 7 

1 \L/3 

l 36 

4 J \ 

. 36 

4 J 


1 - 


1 + 


rr 3-JS7 7 jl/6 ^ 7 3 -JS7 ^76^' 


#17: [s = -2.1795, s = -0.410245 - 1. 74454 - l, s = -0.410245 + 1.74454V] 

These roots indicate that q can not be factorised into linear and quadratic factors with rational coeffi- 
cients. 
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So to find the inverse Laplace Transform: 


ILTC-2.1795) + ILT (-0.410245 - 1. 74454- i.) + ILT(-0 . 410245 + 1. 74454- i.) 

- 2.1795 -t - 0.410245 -t 

2.50551-e - e - (1 . 50551 -C0S(1 . 74454 - 1) + 1 . 23634 - SIN(1 . 74454 - 1)) 


Once more I called the new user contributed function 


2 


\ 


INVLaplace 


3 




s - 4*5 + 2 


2 

+3*s + 5*5 +7 j 


- 0. 4102454876* t 

e • (6*COS(1.744543250*t) + 4*SIN(1.744543250*t)) 


I interrupted simplifying after more than 5 minutes and then achieved a result by approximating the 
expression - but it was quite different to the result obtained by using J V Royle's method (compare the 
two expressions!). I wrote a mail to Terence and asked for advice. His answer was an update of 
LaplaceTransform.dfw from 9 November 2004 which not only simplifies the expression within 
0.77 sec to a bulky expression but also approximates to the same result as given above: 


InvLaplace 


2 

s - 4-s + 2 


2 

+ 3 ■ s + 5 ■ s + 7 j 


- 2. 1795 09024 -t - 0 . 4102454876 ■ t 

2.505504641-e - e ■ (1 . 505504641 -C0S(1 . 744543250 ■ t) + 1 . 256344551 ■ SIN(1 . 744543250 ■ t)) 

- 2. 1795 09024- t - 0 . 4102454S76 ■ t 

Lap lace (2. 5 05 5 04641- e - e ■ (1 . 505504641 -C0S(1 . 74454325 ■ t) + 


1 . 256344531 ■ SIN(1 . 74454325 ■ t))) 

27 2 28 27 

0.002724795640- (3.628195461-10 -s -1.451278183-10 ■ s + 7 . 256390905 -10 ) 

6 7 IS 2 IS IS 

(6.177695-10 ■ s + 1.3464342-10 )■ (1 . 600287999 -10 -s + 1.313021861-10 -s + 5.139696996-10 ) 

(Terence says: Note that the inversion has not returned the original expression. However a plot of the 
original function and the inversion of its Laplace transform show that they are indeed the same func- 
tion. This latest version of LaplaceTransform.dfw is among the files associated to revised DNL#4. 
Many thanks Terence for your immediate reaction.) 


Further extensions of this method are possible: 

(i) Find the formula for the coefficients, A and 5, in terms of P(s ) and Q(s) and their derivatives, 

A B 

for the Partial Fractions 1 — , corresponding to a double root s = a of Q(s) = 0. 

s ~ a (, s-a ) 

Then define a function i lt2 (a) :=A exp(a t) + B t exp(a t) in DERIVE using the for- 
mulae found for A and B. 

(ii) For Laplace Transforms which are not rational functions, finding a Pade approximation valid 
for large 5 and inverting this to get an approximation to the inverse Laplace Transform valid 
for small t. 
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Joint Venture Derive und Cabri 
im Mathematikunterricht 


Herbert Appel und Richard Brand, Schweinfurt 


1. Einleitung 

Das Erscheinen von Programmen fur interaktive Geometrie und symbolische Algebra, die eine neue 
Epoche fur den Mathematikunterricht einzuleiten scheinen, erzeugte in letzter Zeit eine Flut sicherlich 
notiger und guter Artikel (Schuhmann, Engel, Weigand, Weth u.a.) in didaktischen Zeitschriften. Bu- 
cher zu diesem Thema sind momentan noch Mangelware, was sich aber sicher demnachst andern wird. 
Ein Mangel haftet jedoch all diesen Betrachtungen an: zum einen beschranken sie sich meist nur auf 
eine Art von Programm, also entweder auf den Einsatz von z.B. DERIVE oder von z.B. Cabri im MU; 
zum anderen werden darin relativ einfache und grundlegende Beispiele behandelt. Die dadurch mogli- 
che Verquickung von Algebra und Geometrie findet nicht statt. 

Mit der folgenden Aufgabe wollen wir speziell dem Lehrer (direkt an der Front) ein Beispiel an die 
Hand geben, wie er DERIVE und Cabri zur Bearbeitung einer etwas anspruchsvolleren Aufgabe kombi- 
niert in seinem Unterricht einsetzen kann. 

Appearance of dynamic geometry and computeralgebra programs which seem to open a new era for math educa- 
tion initiated a lot of interesting papers in the last time (Schuhmann, Engel, Weigand, Weth ao). All of these arti- 
cles are mostly restricted to only one of these programs and tackle very simple examples. There is very little co- 
operation between algebra and geometry. 

By presenting the following problem we will offer an example how to use Derive and Cabri treating a more de- 
manding task in classroom. 


2. Methodische Hinweise 

Folgende Aufgabe konnte Inhalt einer Ubungsstunde in der 9.Klasse der Sekundarstufe I sein. Sie ist 
beispielhaft fur das induktiv-deduktive Vorgehen beim Finden einer Losung. 

Voraussetzung ist, dass ... 

... die Schuler bereits im Umgang mit der Bedienung der beiden oben erwahnten Programme versiert 
sind oder 

... eine Systemkombination IBM-kompatibler AT + Panel an der Schule existiert. 

Um ein Arbeiten mit Bleistift und Papier kommt der Schuler allerdings schon der Anschaulichkeit 
wegen nicht herum. 


3. Die Aufgabenstellung 


Mit den VektorenH5 und AD n ist eine Schar gleichschenkeliger Trapeze gegeben. 

( 4z-4 ^ 

Dabei gilt: AB = CD n : A( 2 1 1); 5(4.5 | 0.5) und AD n = 

l v 2z + 2.5 y 

Vectors AB and C n D n and define a family of isosceles trapeziums: AB = C n D n : A(2|1); B(4.5|0.5) and 
( 4z-4 ^ 

AD n - 

\2z + 2 .5y 

Fragen / Questions: 


a) Was lasst sich iiber den geometrischen Ort der Punkte D n aussagen? 

b) Welche Kurve beschreiben die Punkte C n l 

a) Which is the locus of points D n ? 

b) Which is the curve described by points C n ? 
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Losungsstrategie: 

Um den geometrischen Ort der Punkte D n zu gewinnen ist eine Parameter-Transformation (mit DERIVE 
oder per Hand) durchzufuhren. Diese ergibt 


AD = 


' x ^ 


x + 9 
\ 2 ) 


; die Punkte D n liegen auf einer Geraden (Frage a). 


Um sich einen Uberblick iiber den geometrischen Ort der Punkte C n zu verschaffen wird mit Cabri das 
(veranderbare) gleichschenkelige Trapez konstruiert. Da in Cabri eine Anfgabe der genauen Punkt- 
koordinaten nicht moglich ist, kann man sich durch Verwendung herkommlicher Zeichenhilfsmittel 
(Zirkel, Lineal und Papier) einen groben Uberblick iiber die Lage der Punkte und Geraden im 
Koordinatensystem verschaffen, eine Idee fur die Konstruktion der Trapeze aus den Angaben holen 
und diese dann mit Cabri durchfuhren. 



Parameter transformation leads to the answer of ques- 
tion a): points D n are lying on a line. 

To obtain an impression about the locus of points C n 
we use Cabri to produce the dynamic variable trape- 
zium. 

We show the “historic” screenshot from 1991 showing 
one of the first German Cabri versions. 


Recent Cabri versions 
give a presentation of 
the locus. We can locate 
pts A and B according to 
the given coordinates. 

You can see that we 
draw the perpendicular 
bisector of AD n . 

D n is moveable on the 
line. 

Reflecting point B with 
respect to this bisector 
gives point C n . 


zu b) Konstruktionsbeschreibung 

Zeichne Gerade g durch A und 5, 

zeichne beliebige zweite Gerade, die g schneidet (Skizze) 

Punkt D n liegt auf der zweiten Geraden, 

zeichne Strecke AD m 

konstruiere Mittelsenkrechte h zu AD m 

spiegle B an h. Es ergibt sich C m 

bewegen des Punktes D n zum Zeichnen der Ortlinie. 
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Man erkennt, dass die Kurve eine unbekannte Form (keine Gerade, Parabel, Hyperbel usw.) ist. Eine 
algebraische Beschreibung der Ortslinie ist sinnvoll. Die gewonnene Ortslinie wird vorsorglich zum 
spateren Vergleich ausgedruckt. 

It is obvious that the curve shows an unknown form (no straight line, parabola, hyperbola, etc). It makes sense to 
look for an algebraic description of the locus, we print the locus for later comparison. 

Zur Veranschaulichung von bekannten und gesuchten GroBen dient die Skizze auf dem Blatt. Die Um- 
formungen werden mit DERIVE durchgefuhrt. Follow the sketch. 



Nach langerer Rechnung ergibt sich der Schnittpunkt der beiden Geraden gnh S(- 31/7; 16/7), deren 
Schnittwinkel a « 37.875° und AS « 6.556 LE. Der Neigungswinkel der Geraden AB ist « -1 1.31°. 


Aus dem Sinussatz folgt: 

6.556 sin37.875° 

o(p) = 

sin(37.875° + P) 

(1) 

auGerdem 

\AB\ = V6T 

(2) 

und 

c(p) = 6(p) + 2 ,45 cos p 

(3) 

und damit C n = (4.5- 

c(P)-cos(11.31° + P), 0.5 + c(P ) • sin(l 1 .3 1° + p )) 

(4) 


We apply the sine rule to express length b depending on angle (3 and find easily an expression for the parallel 
side c((3). Using the coordinates of point B we derive the coordinates of C n which describes the locus in its pa- 
rameter representation with parameter (3. 

To keep work as easy as possible we use the Substitute command three times, which leads us to a huge and 
bulky expression. Students are very happy that they didn't need find it by manual calculation. Derive allows arbi- 
trary long function expressions and is able to plot graphs for functions given in parameter form. 

Um den Schreibaufwand gering zu halten, kommt uns DERIVE mit dem Befehl Substi tute sehr ent- 
gegen. Dreifaches Anwenden ersetzt (1) und (2) in (3) und (3) in (4). 

Heraus kommt ein Riese von einem Term, von dem die Schuler froh sind, dass sie ihn nicht durch 
schriftliches Umformen herleiten mussten (siehe #34 im Listing). 

DERIVE bietet uns nun gegeniiber anderen bisher eingesetzten Kurvenplotprogrammen noch zwei Be- 
sonderheiten: 

- es erlaubt die Eingabe beliebig langer Funktionsterme 

- es kann auch Graphen von Funktionen zeichnen, die in Parameterform gegeben sind. Das nutzen wir 
schamlos aus. 
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Der Vergleich der von Cabri gezeichneten 
Punktmenge und dem Graphen, den uns DERIVE 
zeichnet. ergibt eine ausreichende Ubereinstim- 
mung und ermoglicht dem Schuler eine tiefere 
Einsicht in die Zusammenhange zwischen Al- 
gebra und Geometrie. 

Comparison of the Cabri- and the Derive results 
shows an acceptable correspondence and might 
lead the student to a deeper insight into the connec- 
tions between algebra and geometry. 


Die DERIVE Ausarbeitung / The Derive session 

1 





**28 : 

p_ := ATAN 

— 








5 , 





' b_ SIN(ot) 


P_ 




#24 = 

SOLUE 

= b_ 

#29: 







. SA SINfa + p) 


D 






3 ' s[2G 

#30 : 

-11 .30993247 



#25: 

b_ = - 









7 ■ COS(p) + 9 - SIN(p) 


9 






#31 : 

XC(p) := — 

c(p) ■ C0S(p 

- p_) 



9-^26 


2 



#26: 

b(P) : 

; = 








7-C0S(p) + 9-SIN(p) 


1 






t*32 : 

ycfp) := — 

— + c(p) ■ SINfp 

- p_) 

**27 : 

C(p) : 

:= b(P) + 2 ■ AB ■ COS(p) 


2 




tt33: [xc(p), yc(p) ] 


**34 : 


3 2 

88 ■ COS(p) + 9 ■ COS(p) - SIN(p) (76SIN(p) + 23) 

2 ■ ( 7 ■ COS(p) + 9 ■ SIN(p) ) 


3 2 2 

1 4 ■ COS(p) + 88 ■ SIN(p) ■ COS(p) + 5 ■ COS(p) ■ ( 1 8 ■ SIN(p) + 5) + 99SIN(p) 


2 ■ ( 7 ■ COS(p) + 9 ■ SIN(P) ) 


4. Schlussbemerkung 

Wir hoffen, dem geneigten Lehrer (Leser) einen Eindruck von der Machtigkeit, die eine Kombination 
der beiden Programme darstellt, vermittelt und ihn zum Weiterdenken in dieser Richtung angeregt zu 
haben. 

Ein Programm, das die Fahigkeit von DERIVE und Cabri in sich vereinigt ist unseres Wissens nach 
noch nicht auf dem Markt, aber es wird dem Vemehmen nach daran gearbeitet. 

Concluding comment: We hope that we were able to give an impression of the power which is represented by a 
combination of these two programs and that we inspired some of you to develop more ideas and concepts in this 
direction. It would be great to have a program which combines the capabilities of both tools. 

By rewriting this article from 1991 using the tools of 2004 I found some nice things. 

Please follow the next page(s)! Josef 
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First of all I used Cabri II to find some other loci varying the position of the underlying basic 
lines. 




I sat in the train on my way to give a two-days course on the use of CAS-Tls and DERIVE 
and wanted to tackle Herbert's problem without using trig functions. The coach was full and I 
could not open my notebook. So I took my Voyage 200 and started calculating. 


The idea was to draw a line parallel to AD n through point B. Then draw a circle with center D n 
and radius AB. The intersection point should give C n . Finally I tried to find a reasonable alge- 
braic form for the locus. 


Up! 1 

P^l^.IciS’cWKrKS.Iolci JE’ 1*1 1 

P^l^.IciS’cWKrKS.Iolci JE’ 1*1 1 

t - (2 ■ x - 7) + 14 ■ x - 67 , 

,y= 4 - (t - 2 ) +ba * e 

t-(2-x-7)+ 14-X-67 
y ‘ 4(1 -2) 

■(x - t ) 2 +|y --77 - 9x2 j - 13x2 circ 
Fit. 2 f ?-u-9l -o o ► 

■ cir-c | base + circ 

9 (2 i-(2-x - 29) + 5(2x + 7))(2x - 9) 

lS-Ct- 2) 2 * 

■ zeros(circ , x) + xl 

f 2 ■ t + 5 10 t 3 - 15-t 2 + 3GG ■ t + 109 1 

J 2-t- + 5 10-t^- 15 t^ + 300 t + 109 l 

1 2 2-(5 i 2 - 2 -t + 65) J 

. . , J, i f 2-t + 5 10 t 3 - 15-t 2 t 

■ right base x = < ~ f ► 

{ L 2 2-(5t 2 -: 

ft + S 5t 3 + 2Gt 2 -23t + S9sl 

1 2 2-(5 i 2 - 2 -t + 65) J 

1 2 2-(5-t 2 - 2 -t + 65) J 


...8>/<2*<5*t A 2-2*t+65>>+vt4<t> 


MAIN RAD AUTD PAR EH/HO 

MAIN RAD AUTD PAR EH/H0 

MAIN RAD AUTD PAR EH/HO 


You can see the most important part. Two facts are of special interest: 


o The result is a rational expression and more interesting 

0 We have two loci, one of them appearing as line with the same slope as the locus of D n \\ 

1 plotted all graphs on the Tl using Parametric Mode, which gave a satisfying result. 



It was a long ride, so started Cabri on the Voyage 200 and reproduced my idea using dy- 
namic geometry in full agreement with Herbert Appel's idea of combining the two poweful 
tools and also applying our highly appreciated “Window-Shuttle” principle. 



On the left you can follow 
how to produce the locus. 
There are two intersection 
points! 

One of them giving the trape- 
zium - the other one giving a 
???? 
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Strange things happen plotting the locus of the trapezium point (see right above): There is a 
flaw in the locus and it changes from the curve to the line. This is caused by a deficiency of 
Cabri, because it does not “keep” the right point by passing the position when both quadran- 
gles coincide (in form of a rightangle). Adding the locus of the second intersection point gives 
the missing parts and now we can see both graphs. 




I saved my Voyage 200 session from the Homescreen into a text file appeltxt and hoped 
to be able to tranfer the calculation to Derive 6! I added some comments in the text file and 
closed with four lines which should give the parameter form in Derive. Plotting them in the 
2D-Plot window should immediately give the graphs. 




■ Comand U i ew Execute 


k r 


C: [[i][U2+9>2]]+d 
C: t+xt2(t) : tx2+9x2+yt2(t) 

: direction vector of II line 
C:d-a 

C:y-lx2=(tx2+7x2)x(t-2)=Kx-9x2) 

C:solve(y-lx2=(t+7)=K2*x-9V(4*(t-2)) ? y> 

C:y=<t*<2*x-7) + 14*x-67V<4*<t-2>)+base 
C: fx-t^2+Cy-tx2-9x2^2-13x2^circ 
C: circ I y=(t=K2*x-7) + 14*x-67)x( 4*(t-2) >+c 
ire 

: Intersect line and circle 


I^UTandl 

f Y FH T FE 1 

|Uiew|Execute Find...| 


C: C5+t^3+20+t^2-23*t+898Vt:2*C 
_65myt4(t> 

: parameter representation! in 

C: Extl(t) 7 ytl(t) 1 
C: Ext2( t) 7 yt2(t) 1 
C: Ext3( t) ? yt3(t) 1 
C: Ext4( t) 7 yt4(t) 1 

5*t^2-2*t+ 

Derive 

MAIN RAD AUTD PAR 


Having arrived in my hotel room in St. Johann/Salzburg I opened the notebook, connected it 
with the V 200 and launched Derive 6.10 and hoped that InterConnectivity would do it. 


I issued File>Ti Handheld>import from and typed in appeltxt. Transfer started and: 


Intersect line and circle 


# 29 : xl := SOLUTIONS 


9 . ( 2 . t- ( 2 -x - 29 ) + 5 . ( 2 -x + 7 )).( 2 -x - 9 ) 5 -t 5 .t.( 2 .x - 5 ) 


16 - (t - 2 ) 


5 -x 


43 .x 


845 

16 


x, Real 


# 30 : 


xl := 


I 2 I 3 2 

18 . 1 1 - 5 ■ t + 6 | + 10- 1 + 3 • t + 210- 1 + 217 


2 

2 - ( 5 • t - 2 * t + 65 ) 


I 2 I 3 2 

18 . 1 1 - 5 • t + 6 | - 10- 1 - 3 • t - 210 . t - 217 


2 

2* (5 - 1 - 2- 1 + 65 ) 


# 31 : yl := RHS 

r 

VECTOR 

r 

base, x, 

3 2 

2 -t + 5 10 - t - 15 - t + 300 - t + 109 

>> 

2 ' 2 





2 ■ ( 5 ■ t - 2 -t + 65 ) 

Jj 


It seems to work. I receive the correct - both loci and should be satisfied (next page). In- 
specting the Derive procedure more accurately you will find quite another result for the solu- 
tions of the important equation which delivers the x-values of the intersection points (see 
expr #30). 
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parameter representation in Derive 
#49: [[xtl(t), ytl(t)]] 


#50: 


#52: 


#54: 


5 ■ t t 

+ 2, 1 

2 2 


#51: [ [xt2(t) , yt2(t)]] 

t 9 

2 2 

#53: [ [xt3(t) , yt3(t)]] 

2-t + 5 t + S 


2 2 
#55: [ [xt4(t) , yt4(t)] ] 



3 2 

10-t - 15. t + 300-t + 109 

3 2 

5-t + 20- t - 23. t + 898 

2 

_L 2- (5-t - 2-t + 65) 

2 

2- (5-t - 2-t + 65) J. 


Removing the abs from expression #30 and simplifying once morel obtain the same expected 
expressions. (See the end of the page for another way to obtain the expected result.) 

2 3 2 2 3 2 ' 

18. (t - 5-t + 6) + 10. t + 3-t + 210- 1 + 217 18- (t - 5-t + 6) - 10. t - 3-t - 210. t - 217 

2 ' 2 
2* (5 • t - 2*t + 65) 2* (5 . t - 2-t + 65) 

3 2 

2-t + 5 10. t - 15. t + 300-t + 109 

2 2 

2.(5.t - 2-t + 65) 

In this case InterConnectivity was very helpful for me. However, it must be said that this tool 
is not perfect at the moment. You have to consider carefully how to use the |-operator on the 
Tl to have a correct translation in Derive. I was said - and from my experience with the de- 
veloppers I have no reason to doubt - that the programmers are busy to improve InterCon- 
nectivity between the CAS-Tls and Derive. 


Having used successfully InterConnectivity I wanted to introduce slider bars to have a dy- 
namic Derive which is very close to Dynamic Geometry and forms a useful link between al- 



This is a way to overcome Derive's solution of equation #29 including the absolute values: 
Factor the expression and then investigate the factors to get the two solutions. 

2 3 2 

(2-x - 2* t - 5). (2-x. (5. t - 2-t + 65) - 10-t + 15. t - 300-t - 109) 


16- (t - 2) 

' 2-t + 5 ' 

S0LUTI0NS(2.x - 2-t - 5 = 0, x) = 


2 3 2 

S0LUTI0NS(2.x ■ (5 • t - 2-t + 65) - 10-t + 15. t - 300-t - 109, x) 


3 2 

10. t - 15. t + 300. t + 109 
2 

2. (5-t - 2-t + 65) 
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ITERATION 


Module on Sequences and Series 
with Derive , Part 1 

David E. Hodgkinson, Liverpool 


Iteration is a method of repeated calculation that has become very important because computers 
can iterate quickly and accurately. The process of counting can be represented as an iterative process 
by starting at 1 and repeatedly adding 1. One mathematical way of writing this is /?:= 1 and n <— n + 1 
This means that n is initially assigned to 1 then repeatedly increased by 1. Usually a way of stopping 
the repetition is included, either by saying how many times the assignment is repeated or providing a 
condition to stop the repetition. 

DERIVE has various functions that provide iteration facilities. To produce a list of the first twelve 
positive integers use the function ITERATES with four arguments. Commas separating the arguments 
of a function must be typed. Hence typing the following instruction should produce the associated list 
of integers. 


#1: 

ITERATES (n + 1, 11 , 

1, 

ID 

#2: 

[1, 2, 3, 4, 5, 6, 

7, 

8, 


9 , 10 , 11 , 12 ] 


is obtained because the second argument of iterates is updated by the assignment to the first argument 
i.e. n <— n + 1, the third argument is the initial value of n and the fourth argument is the number of 
iterations required. Hence changing 1 to 0 in iterates gives a list of 13 consecutive integers beginning 
with 0. It is fairly obvious that as the number of iterations increases the final value in the list will in- 
crease, this is an example of a divergent iteration. 


Another example is the assigment 



the first six values of this iteration are found from 




r i 

2 ' 


#3: 

ITERATES 

^ 2 

x + — 
x ^ 

, x, 1, 6 

J 


3 17 577 665857 886731088897 1572584048032918633353217 

2 12 408 470832 ’ 627013566048 ’ 1111984844349868137938112 


However, the list of longer and longer fractions is not very informative. Since iteration is mainly used 
for numerical calculations, to get decimal approximations use the approx-button (»): 

A5: [1, 1.5, 1.416666666, 1.414215686, 1.414213562, 1.414213562, 1.414213562] 

i.e. approximating the fractions (or see below for how to repeat the procedure in approximate mode) 
shows that the iteration is settling to 1.414213. This is recognisable as V2 to 6 decimal places. Notice 
that this iteration is approaching a limiting value as the number of iterations increases, this is an exam- 
ple of a convergent iteration. 
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2. Try the iteration x n+l = — 


x +- 


N 


for TV = 27, 8,4, 2, 3. (x 0 =l). 


x n J 


What value does this iteration converge to for given TV? 


3. Try the iteration x = 2x -TV-x 2 

J n + 1 n n 

where TV = 2, 3, 4, 5, 6 and the initial value of x lies between 0 and 1 . What is the value of this it- 
eration for different TV? 


4. Find V2 for 10 and 20 decimal places to see the speed with which the convergent answer is 
reached. If you have been in approximate mode, do not forget to return to exact mode using the 
above instructions. 

What value does this iteration converge to for given TV? 


SEQUENCES AND SERIES 

In the previous examples one list had the values 1,2, 3, 4, whilst the other had to 4 decimal places 
the values 1, 1.5, 1.4166, 1.4142, 1.4142. Lists of numbers separated by commas are called sequences. 
Mathematically speaking sequences are functions with domain the set of integers, but instead of using 
the standard/^) or g(n) notation, the alternative x n or u n are preferred, but with DERIVE the standard 
notation will be used. 

Closely connected with sequences are series. Series are obtained by adding together the terms of a 
sequence. The value of this combined addition is called the sum of the series. 

e.g. the finite sequence 1,2,3,4,5,6,7,8,9,10, 

gives the series 1+2 + 3+ 4 + 5 + 6 + 7 + 8 + 9+10, 

and the sum af the series is 55. 


Using Derive the sum can be obtained iteratively by using an array v in which the first element is 
the successive numbers of the series and the second element is the sum of the series. The iterative 
formulae being used are x:=l and s:= 1 initially and x is increased by 1 each time whilst s is in- 
creased by x i.e. x^x + 1 , s <— s + x. 

The initial value of v is [1,1]. The first element of v represents the next number in the series, the 
second element represents the current value of the sum. Derive updates its copy of v using the function 
ELEMENT, which extracts elements of arrays. The next number in the series is found by using 

ELEMENT (v , 1)+1 

Works still but more comfortable use v sub 1 + 1 which appears as ^ + 1 
The value of the sum is found from 

l+element(v,l)+element(v, 2) orbetter 1+v sub 1+v sub 2 = 1+v +v . 

1 2 

The first two elements represent the new value of the nth term added to the third term representing the 
previous sum. 
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The command to automate this is: 

ITERATES (T 1 + v , 1 + v + v 1, v, [1, 1] p 6) 

#19 : L 1 1 2 J 

r 1 11 

The answer is a two column array in which the first contains the terms 
of the series and the second column the corresponding sum of the 
terms so far. 

5 15 

6 21 

7 28 

Naturally infinite sequences and series do occur. The idea of convergence can be used with series 
as well. If the sum to infinity is finite then a series is called convergent . The problems of investigating 
the convergence of series are many and varied and still the subject of active research. This module will 
give some ideas of the use and problems with sequences and series. 

First note well that a convergent sequence does not mean that the corresponding series will con- 
verge . 

The infinite sequence 1, 1, 1, 1, 1, 1, 1, 

is convergent, but the infinite series 1 + 1 + 1 + 1 + 1 + 

is divergent. 


2 3 

3 6 

#20: 4 10 


Do not be deceived by the size of the first few terms in a series, it is the long term behaviour of a 
series that matters. 

Use Derive to find the sum of the first few terms of the series 


100 ! 100 ! 100 ! 100 ! 100 ! 100 ! 



1205 08715 72087401733412468141 767066166 3407403071516651940087400101025 5 5 5965 85 06311682 76210090701 0466871597242370877488r- 
54379741889251749 B964 80000000000000000000 

The sum looks very large and impressive. However, this series has a finite sum, since the denominator 
eventually gets much larger than the numerator. Try 

100 ! 

100 

100 

100 ! 

200 

200 

100 ! 


300 


300 


and simplify the expressions. 
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Whereas, the sum of the series in which each term is the reciprocal of the corresponding one in 
this series, starts small, but eventually grows without bound. Repeat the trial by looking at the first few 


sums of 


100 ! 


Then examine the later terms, these are very large and increasing so that the sum would 


grow without bounds i.e. the series is divergent. 


However, you must not be mislead into thinking that individual terms must be large for the series 
to be divergent, it is the cumulative effect of adding an unlimited number of terms that is important. In 
fact the series 


1 + — + - + — + - + - + — + 

2 3 4 5 6 7 

caused by adding the reciprocals of the integers is divergent, even though the individual terms get very 
small. It is impossible to show the divergence by just finding the sums to a finite number of terms, the 
divergence needs a theoretical proof. This series is a famous one in mathematics and is called the 

harmonic series. 


ARITHMETIC SERIES 


The series 

1 + 2 + 3 + 4 + 5 + ... 

is obtained by summing the terms of the sequence 

u n +\ = u n + 1 with it i = l. 

In fact any sequence of the form 

u n +\ = u n + d 

with u\ given an initial is called an arithmetic sequence and the sums of the terms of these sequence 
are called arithmetic series. 

One of the fascinations of mathematics is finding alternative ways to solve problems. Faced with the 
problem of adding 

1+2+3+4+5+6+7+8+9+10 

the obvious way is to do the addition and get 55 either with pencil and paper or using the iteration and 
calling 9 times. Summing from 1 to 100 is much more arduous and all programming languages have 
ways of instructing the computer to carry out repetive steps. An alternative way is to find a method of 
getting the answer from a formula and avoiding the repetetive steps. This is usually not possible, but in 
most instructive examples the alternative of finding a formula is possible. 

In the case of arithmetic series, this alternative of a formula for the sum is possible. The following 
trials are aimed at producing a formula for the sum of an arithmetic series. 

Considering the sequence u n+ \ = u n + 1 with u x = 1 

it is very easy to find the 10th term, it is 10, the 100th term is 100 etc. so in this case the nt\\ term is n. 

The sum of the terms is less obvious. So try to author a few iterations consulting the previous ex- 
ample if you need to see how. The answer gives two columns: in the first one are the values of n and 
in the second one are the values of the sum to n terms. So for the n sequence 

1,2, 3, 4, 5, 6, 7 


the s (i.e. sum) sequence is: 1, 3, 6, 10, 15, 21, 28. 
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To try and work out a formula connecting the sum to n terms with the value of n , trial and error 
can be tried, but a systematic approach is best. Assuming s = an + b i.e. a linear relationship, fit this 
equation to the first two terms to find values of a and b see if the formula predicts the third and higher 
values. When n = l, s= l;n = 2,s = 3. DERIVE expects you to enter linear equations in the form of a 
vector, so the [,] brackets must be typed, around the equations. 


[I = a + b p 3 = 2-a + b] 

S0LVE([1 = a + b p 3 = 2-a + b] p [a, b]) 
[a = 2 a b = -1] 


Hence author s = 2n - 1 

To test this conjecture that this is formula is true if n = 3 gives s = 6. Substitute for n and 5 : 

s = 2 mi - 1 The linear guess is wrong. (Undoubtely you will not need a CAS for this calcula- 
tion, but it might be useful for later to learn this procedure.) 


6=5 

Trying again with a second degree equation s = an 2 + bn + c there are three unknown coefficients a , b 
and c, so the first three terms will be needed and then the predictions for the fourth and higher terms 
can be attempted. Author the first three equations as a vector and solve, remember to type the [ ] 
around the equations. 

SOLVE([a + b + c = 1, 4 a + 2 b + c = 3, 9-a + 3-b + c = 6] p [a, b P c]) 

1 1 

a = a b = a c = 0 

2 2 

n 2 n 

and these answers make 5 = — + — . To test for higher terms it is possible to use a DERIVE command 

which allows you to do repeated substitution into a formula. It is called VECTOR and has four argu- 
ments, the first for the expression, the second for the variable, the third for the starting value and the 
fourth for the finishing value. 


( 2 


VECTOR 


11 


k 2 



Ip 


\ 


7 = [1, 3, 6 p IQ, 15 p 21 P 28] 


giving complete agreement. 


Exercises . 

1 . Consider the iterative sequence 

u n + 1 = u n + 2 with u i = l. 

Generate the first seven terms of this sequence and also the sum for n = 1, 2, 3, 4, 5, 6, 7. 
Use your results to produce a formula for the nth term and a formula for the sum to n terms. 

2. Use the iterative formula Uyi -\- 1 Mfi ~ E d • 

Choose your own u x and d and see if you can find the formulae for the nth term and the 
sum to n terms. 
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Working with the general arithmetic sequence u n+ \ = u n + d,u\= a , you can find its sequence of 
partial sums by working with iterates and a vector with three elements. This is so that the first element 
can keep count of the number of terms, the second element can update the current value of the nth 
term and the third element can add the nth term to the previous sum. 

ITERATESCT1 + v , v + d, v + v + d], v, [1, a, a], 7) 

L 1 2 3 2 J 


The first column sequence represents the 
number of terms, the second column sequence 
gives the value of the nth term ( u ) and the third 
column the sum sequence gives the value of the 
sum to n terms ( 5 ). The problem is to find a for- 
mula relating the value of n with the values of the 
nth term u and the sum 5 . 


la a 

2 a + d 2-a + d 


3 

a 

+ 

2 ■ d 

3-a + 3-d 

4 

a 

+ 

3-d 

4-a + 6-d 

5 

a 

+ 

4 ■ d 

5 ■ a + 10 ■ d 

6 

a 

+ 

5-d 

6 ■ a + 15 ■ d 

7 

a 

+ 

6 ■ d 

7- a + 21. d 

S 

a 

+ 

7 ■ d 

S-a + 2S-d 


Can you spot the relationship between n and ut 

In conventional work the first term is usually denoted by a and the commom difference by d, so 
that the formula for the nth term would be a + (n-l)d, which should be equivalent to the one you 
obtained. 

The connection between n and s is not as obvious, since both a and d have changing coefficients. 
Can you spot the formula connecting the nth term with the coefficient of al Write down the sequence 
of coefficients of d. Have you seen this sequence before? If so can you construct the formula relating 
the coefficient of d with nl 

Finally can you combine your two results to show that your formula is eqi valent to 

s n = ~' \^ a + (n - 1) • d\ 

(to be continued) 


(David E. Hodgkinson, Department of Applied Mathematics and Theoretical Physics, University of Liverpool) 
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How to write my own DEMO-file? 

Josef Bohm, Wiirmla 

When I want to use DERIVE for demonstrating the procedure how to solve a problem I ? m glad to have 
of a Demo-file at my disposal. Students or colleagues start with File>Load>Demo File>Filename. 
They are able to stop the ’’demonstration” at any position to do their own calculations, change the op- 
tion settings, plot partial results etc. In addition they can find in the message line a short description, a 
little hint, an important explanation, 

I was curious and inspected an original DERIVE -.dmo-file. In DOS: type arith . dmo. It's pure 
ASCII-Code. I saw the message to be printed on screen in the first line, starting with a semicolon. The 
next line was a DERIVE -expression, written in its very own syntax. That was all! 

Now let me write a Demo-file called SKETCH . DMO to demonstrate investigating a function graph. 
You can use each editor, which creates ASCII-files. I prefer WORD and save my file as ASCII-file 
with the extension .dmo. Another way to edit this file is in DOS: 

copy con: sketch, dmo etc. 

;This is the function to be investigated 
f (x) :=(x A 4-llx A 3+29x A 2+35x-150) /30 
;f(x) = 0 to obtain the zeros 
SOLVE ( f (x) =0 , x) 

;f ! (x) is the 1st derivative which is defined as fl (x) 
fl (x) :=DIF (f (x) ,x) 

; zeros of f 1 (x) give the relative extrema, ESC, Approx, Notation 
Decimal, Simplify 
SOLVE (fl (x) =0,x) 

; These are the turning points! Set back to Exact and Rational 
[ 5 , f (5) ; -0. 47038, f (-0.47038) ; 3. 72038, f (3.72038) ] 

; Hurry up for the 2nd derivative 
f 2 (x) : =f 1 1 (x) 

;f IT (x) = 0 gives the inflection point(s), ApproX, Notation Decimal 
SOLUTIONS ( f 2 (x) =0, x) 

;Here are the complete coordinates! Study the Derive command! 

VECTOR ( [ u , f ( u ) ] ,u, SOLUTIONS ( f 2 (x)=0,x) ) 

;We need the slope in the inflection points for a fine sketch 
[fl (1.098) , fl (4.402) ] 

; Let's proceed with the indefinite integral in Exact Precision 
INT (f (x) , x ) 

;Whatis the area between graph and x-axis between the zeros 
ABS (INT (f (x) ,x, -2, 3) ) +ABS (INT ( f (x) , x, 3 , 5 ) ) 

;ESC, Change to 2D-Plot Window, highlight f (x) and plot the graph 
f (x) 

/Finally we will shade the area 

( -2<=x<=3 AND f(x)<=y<=0) OR (3<=x<=5 AND 0<=y<=f (x) ) 

This was from 1991. Nobody will work with copy con .... in 2004?? 

We have other choices: 

(1) If you like the nostalgic way you might work with the DOS-Editor. So switch to DOS and type 

edit sketch. dmo. 

(2) You can do it from the DERIVE 6 environment or 

(3) You can use your handheld CAS-Device and apply InterConnectivity. 


D 


-N-L#4 


Josef Bohm: How to write your own DEMO-file! 


p35 


Don't use the textprcessor now because of differences between ASCII and ANSI-code! 

This is the first - nostalgic - way. You must be very careful typing the DERIVE commands in the cor- 
rect syntax and typing them in uppercase. In DOS mode call edit sketch . dmo and start editing: 



;This is the function to be investigated 
f(x):=(xA4-ll*xA3+29*xA2+35*x-150)/30 
;f(x) = 0 to obtain the zeros 
solve(f(x) , x) 

;f‘(x) is the 1st derivative which is defined as fl(x) 
fl(x) : =DIF(F(x) , x) 


Having finished and saved the file you should be able to load this file as a demo file in Derive. The 
first expression appears in the Algebra Window and in the message line on the bottom you will find 
your comment. 




1^31' El 

B S U & t t J 


#1: f(x) 


4 3 2 

x - 11- x + 29 -x + 35-x - ISO 


30 


The message line (bottom left): 


# 2 : 


4 3 2 

x - 11- x + 29-x + 35-x - 150 

f(x) := 

30 


This is the function to be investigated 

Start | §T] dnl04rev - Microsoft Word | Paint Shop Pro - [Image... 


You can write your complete DEMO-file with the Editor. But possibly you will prefer the DERIVE 
environment, because you want to convert a successful DERIVE-SQSSion into a Demo file and without 
typing in all the stuff once more in the editor. 

So do the following: 


Highlight the expression which should be executed, go for Edit>Annotation and write your com- 
ment as shown below. Later this comment will appear in the message line. 


#9: 


# 10 : 


5 

f<5 

-0.47038 

f (-0.4 

. 3.72038 

f (3 . 72 

5 


-0.47038 

-5.295: 

. 3.72038 

0.2249: 


Annotation for expression ttl 1 


El 


Text 


Hurry up for the 2nd derivative! 


OK 


Cancel 


# 11 : 


I recommed to delete all simplified or approximated results before saving the file. So I would delete 
expression #10. Otherwise they would appear twice. 
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When you are ready with the session then you have to save the file in such a way that DERIVE will 
recognize it as a demo-file. 

You must take care of two details: 

• Save the file as a Math- (mth-) file 

• Type the name of the file under quotes with the filename extension . dmo. 



If you are the happy owner of a TI-92 PLUS, Voyage 200, TI-89 or TITANIUM then you have a third 
chance to produce a DEMO-file: 


At first create your demo as a script on the hand- 
held device. 




■ Command U i ew Execute 




:This is the function to be investigate 

C: < x A 4-l lx A 3+29x A 2+35x-15G V3G+f 
:ftx>=G to obtain the zeros 
C:soluetftx>=G ? x> 

is the 1st derioatioe which is d 
efined as flCx> 

C:d<f<x>,x>+fKx> 

:f'(x>=Q leads to the extremal oalues, 
ESC, Approx, Decimal, Simplify 
C:soluetfKx>=G,x> 

: These are the extremal oalues, Notatio 


Import from Handheld 


Data from 

Tent Variable (File): (sketch 



Send to 

E Worksheet (Automatic simplification) 
E Worksheet [No simplification) 

F DMO File |sketch| 



It is not necessary to save the *.dmo-fIle under the 
name of the script! ! 


OK 


Cancel 


All your one-line comments in the script which are before an expression will appear as comments in 
the DEMO-file. 


If you would import the same file as a worksheet 
then you would find your comments in text- 
boxes. 

The picture shows the start of the file. 


This is the function to be investigated 


#1: f(x) 


4 3 2 

x - 11-x + 29-x + 35-x - 150 


30 


f(x)=0 to obtain the zeros 

#2: SOLVE (f 00 = O, x, Real) 
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CALL FOR APPLICATIONS 

International Spring School on the Didactics of Computer Algebra 

April 27-30, 1992 
Krems, Austria 

Computer Algebra is a new technology with an i mm ense impact on techniccal sciences. The 
personal computer algebra system DERIVE has found its way to all levels of education. It has 
pleased many teachers but also surprised others. The availability of DERIVE on very inex- 
pensive hardware and the new palmtop size Pcs enables more and more students to use the 
system for their assignments or in school. 

How will math teaching need to be changed to cope with thie new development? this spring 
school aims to bring together people with an interest in developing a didactics of computer 
algebra for the math education of students aged 12-18. 

The program will consist of presentations of classroom experiments, discussions and working 
groups. The organizers provide full lodging and boarding. The number of participants is 
strictly limited. 

We invite your application to participate. Applications should include the following items: 
full name and address, 

list of activities related to 'computer agebra in education' (own experiments, 

publications, workshops, seminars, etc.), 

an extended abstract on what you will present at the meeting, 

to be sent to the following address: 

Research Institute for Symbolic Computation 
Univ. of Linz, att. Dr. B. Kutzler 
A-4040 Linz, Austria, fax +43-7236-3231-30 

to arrive no later than 

February 15, 1992. 

The Notification of acceptance will be mailed February 29,1992. 

This spring school is jointly organized by the following institutions: 

Landesschulrat fur Niederosterreich, Wipplingerstr. 28, 

A-1010 Wien, Austria (LSI Dr. Helmut Heugl) 

¥ Bundesministerium fur Unterricht und Kunst, Wipplingerstr. 20, 

A-1010 Wien, Austria (MR Dr. Eduard Szirucsek) 

¥ Research Institute for Symbolic Computation, University of Linz, 

A-4040 Linz, Austria (Dr. Bernhard Kutzler) 





D 


-N-L#4 


More Graphics from Mr Serif's Treasure Box 


p38 





